Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 20 January 2013 (MN IATeX style file v2.2) 



Fluid-like entropy and equilibrium statistical mechanics of 
self-gravitating systems 



Dong-Biao Kang* and Ping He 



Oh 

< 

m 
(N 

o 

u 

43 
6 

(N 

> 
in 

On 

m 

O 



is 



STey Laboratory of Frontiers in Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Science, Beijing 100190, China 



20 January 2013 



ABSTRACT 

The statistical mechanics of self-gravitating systems has not been well understood, and still 
remains an open question so far. In a previous study by Kang & He, we showed that the 
fluid approximation may give a clue to further investigate this problem. In fact, there are 
indeed many dynamical similarities between self-gravitating and fluid systems. Based on a 
fluid-like entropy, that work explained successfully the outer density profiles of dark matter 
halos, but there left some drawbacks with the calculation concerning extremizing process 
of the entropy. In the current paper, with the improved extremizing calculation - including an 
additional differential constraint of dynamical equilibrium and without any other assumptions, 
we confirm that statistical-mechanical methods can give a density profile with finite mass and 
finite energy. Moreover, this density profile is also consistent with the observational surface 
brightness of the elliptical galaxy NGC 3379. In our methods, the density profile is derived 
from the equation of state, which is obtained from entropy principle but does not correspond 
to the maximum entropy of the system. Finally, we suggest an alternative entropy form, a 
hybrid of Boltzmann-Gibbs and Tsallis entropy, whose global maximum may give rise to this 
equation of state. 
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1 INTRODUCTION 

We know that two-body relaxation may play a crucial role in sys- 
tems such as stellar clusters and the center of galaxies, while plays 
no role in collisionless systems such as galaxies and dark mater ha- 
los, whose time-scale s of two-body rel axation are longer than the 
age of the universe. iLvnden-Belll d 19671) proposed 'violent relax- 
ation', which is another relaxation mechanism due to the rapidly 
changing gravitational field when systems collapse. Violent relax- 
ation can allow the collisionless system to be in equilibrium. If 
we assume that the Boltzmann-Gibbs (BG) entropy is applicable 
to self-gravitating systems 



/In/dr, 



(1) 



where / = f(x, v) is the distribution function, and dr = d 3 xd 3 v is 
the phase-space element. By maximizing the BG entropy with the 
constraints of given mass and energy, we can obtain the Maxwellian 
distribution for the system, which corresponds to an isothermal 
density profile and hence cannot be accepted because of its infi- 
nite mass and infini te energy. Then, to eliminate this inconsistency, 
lLvnden-Belil ( ll967l) also proposed the incomplete relaxation, which 
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means that the isothermal solution is only valid for the central re- 
gion. 

In fact, t he rea son for this inconsistency is a result pointed out 
bv lAntonovl ( ll962h that the maximum of the BG entropy of self- 
gravitating systems does not exist, so t he thermodynamic equilib- 
rium state may not exist. Binney (see iBinnev & Tremainell2008l) 
also proved that if one can divide the system into a central main 
body and an outer envelope, then the BG entropy of self-gravitating 
systems can increase without upper bound by increasing the sys- 
tem's central concentration. So there is no global maximum of en- 
tropy for self-gravitating systems, or we can question the validity 
of the BG entropy for self-gravitating systems. 

Despite these difficulties, statistical-mechanical methods have 
always been explored in an attempt to expla in the properties of self- 
gravitating systems. Trem aine et alj {l986) discussed the general- 
ized H -function 



H[f] = 



C(/)dr, 



(2) 



where C(f) is any convex functional of /. In fact, the BG entropy 
(Q3 is just one of such C(/)s and may be only valid for short- 
range interacting systems. F or long-rang e interacting systems, such 
as self-gravitating systems, iTsallisI dl988l) generalized the form of 
BG entropy in mathematics as 
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s q [f] 



(f - f)dr, 



(3) 



and established the non-ext ensive statistica l mechanics. When q = 
1, it recovers to BG entropy. Hansen (2004) used this entropy to in- 
vestigate the density slope of dark matter halos. However, the max- 
imum of this entropy cannot give a correct density profile of self- 
gravitating system. Some authors also studied the ther modynam- 
ics an d kin etic equat i ons of long-range systems, such as lChavanisl 
d2006h and lChavanij (l2010h . 

Besides, some works also applie d some sp e cial en tropy forms 
to self-gravitating systems. Iwhite & Naravanl ([l987) (hereafter, 
WN87) considered the specific entropy form of ideal gas. With 
fixed mass, fixed energy and finite extent, WN87 seeks the equilib- 
rium configuration of elliptical galaxies by maximizing the entropy 
of ideal gas, but the result seems to be not good. However, in this 
paper we wi ll argue abou t the validity of their calculatio n. Our pre- 
vious works lHe &Kand feOld) and llCang & Hel 1201 ll ) (hereafter, 
KH11) employed this entropy form to study the equilibrium con- 
figurations of dark matter halos. Moreover, KH11 considered the 
general anisotropic situation and found an equation of state whose 
density profile has finite mass and extent. Here we will also discuss 
main properties of this density profile. 

In this paper, we demonstrate that the results of WN87 and 
KH11 should have no essential differences and propose another 
kind of entropy form including both BG entropy and Tsallis En- 
tropy. This paper is organized as follows. In section [2] we em- 
phasize the dynamical similarities between self-gravitating systems 
and fluids. We restate the work of KH1 1 with some differences and 
mainly compare the results of WN87 and KH11. In section [3] we 
propose the new entropy form of the hybrid BG and Tsallis entropy. 
We present the conclusion of our work in section[4] The appendixlAl 
explains the derivation of the new equation of state which appears 
in sectionl2?2l 



2 BASIC THEORY AND FORMULAE 

2.1 Connections between fluids and self-gravitating systems 

For collisional system s , some authors (e.g., iLarsonl Il970t 
iLvnden-Bell & Eggletor] Il980l : iBettwiesej 1 1983b integrate the 
Boltzmann equation to obtain a set of moment equations, which 
can be solved to follow the evolution of the system. Under some 
reasonable assumptions, only the low-order moment equations are 
reserved and become fluid-like conservation equations. So the col- 
lisional gravitating systems may be approximated by fluids whose 
hydrodynamic equilibrium can be shown by 



(4) 



where $ and p are the gravitational potential and pressure respec- 
tively, and the density p satisfies the Poisson equation 



V 2 $ = 47rGp. 



(5) 



While the collisionless system is described by the phase-space 
density / = f(x,v) wh ich satisfies the collision less Boltz- 
mann equation (CBE, see iBinnev & Tremainel 120081) . For sim- 
plicity, we restrict ourselv e s to the spheri cal systems in this 
work. ThenlAhn & Shapiro! d2005l) (also see iTevssier et al.lll997l ; 
Subramanian et al .120001) proved that with the three assumptions of 
spherical symmetry, skew-free velocity distribution and isotropic 
velocity dispersion, we can also integrate CBE as the case for col- 
lisional systems, and the derived moment equations are the same 



as fluid conservation equations for a gas with adiabatic index 5/3. 
When the collisionless systems are in dynamical equilibrium, we 
can get the Jeans equation from CBE: 



d(pa 2 ) 



■2/3 



(pa r 



= -P- 



dr r dr 

where p satisfies equation ((5J, with 



p(x)= / /(*,v)d J v, 



(6) 



(7) 



and /3(r) — 1 — of /2<r 2 being the anisotropy parameter, and of and 
cr 2 are the tangential and radial velocity dispersion, respectively. As 
in KH1 1, we define the effective pressure as 



Then equation (|6j can be rewritten as 



dP 



GM 



(8) 



(9) 



dr ^ dr / - 

which is formally equivalent to equation (|4}. Hence, with the def- 
inition of the effective pressure P, such a formal equivalence sug- 
gests that there should be some dynamical similarities between self- 
gravitating systems and fluids. 



2.2 Equation of state 

We will reproduce the equation of state in KH1 1 with the improved 
methods of the variational calculus. Based on the analysis of the 
above section, we infer that the fluid-like entropy with the adiabatic 
index 5/3 is really applicable to current self-gravitating systems as 



ft 



ps47rr dr : 



p3/2 

pln( — r r—)4Tvr dr, 



(10) 



where s = ln(P 3/2 /p 5/2 ) is defined in KH11. Setting y(r) = 
M(r) = J Q r 47rpr 2 dr, the mass function of the system, we obtain 



P = 



V 



4 7rr 2 



(11) 



The mass constraint for extremizing the entropy is equivalent to a 
fixed end-point condition as 



M t = y{r)\ r=00 . 
The energy constraint is 

E t = E k + E v = i / po 2 tot A\ - 



(12) 



2 r — r'\ 



... 1- / ( (3 W) pair 2 -GMpr)Ar 



(13) 



in which we used a change of integration order to derive the first 
term in the last line of the above equation. It is different from KH1 1, 
in that (1) we have not assumed that the system is in hydrostatic 
equilibrium, so the virial theorem cannot be used; (2) P is treated 
as an independent variable rather than a known function. Then the 
constrained entropy will be 



St,i — St + \Et 



(14) 



where A is a Lagrangian multiplier. Substituting equations | |10) , 
i fTTT l and fi"3l into dT4t . we obtain 
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St,r — 



/ F(r,y,y',P,P')dr 
Jo 



(15) 



We use entropy principle, i.e. SSt, r = 0, to extremize St, r and 
obtain the Euler-Lagrangian equation for P and y respectively, as 



P = 
and 



-AP, 



(16) 



-A 



(17) 



dln(P*/pi) 
dr 

In combination with equation {5), the solution of equations l |16t and 
l |17t is an isothermal sphere, which does not require that j3 = 0. 
However, this solution cannot be accepted, because (1) as men- 
tioned previously, it has infinite mass and energy; (2) it has only 
one shape parameter A but there are two constraints, and thus this 
solution is incomplete. Nevertheless, if we discard equation d!6t 
and only substitute equation (|9) into the RHS of equation (T7J, we 
can get the result of KH1 1 



-AP + vP 



.3/5 



(18) 



where v is an integral constant. 

Here we summarize the main properties of equation ( 118) 
shown by KH1 1: 

(i) There are two shape parameters A and v that can be adjusted 
to make equation {\%\ satisfy the mass and energy constraints, and 
the existence of v term amounts to introducing a 'truncated' radius, 
which ensures the finite mass of the gravitating system. 

(ii) A must be negative and the v term may play an important 
role only at large radii. The density profile is shown in Figure Q] 
which is similar to a 'truncated' cored isothermal sphere. We can 
get different core radii and the 'truncated' radii by specifying dif- 
ferent values of A and v. As shown in Figure [T] this density pro- 
file is only consisten t with the dark matter halo's Einasto profile 
dNavarro et aTll2010t) in the range of r/r_2 > 0.2. Nonetheless, 
this solution should be acceptable. 

According to these results, we can identify an interesting 
structure, i.e. index-varying polytropic system, with the polytropic 
index (see equation J28t and below) changing from 1 to 5/3 with r 
increasing from to oo, and correspondingly, the state of the sys- 
tem varies from isothermal to adiabatic, which is different from the 
'core-halo' structure. However, this structure may not correspond 
to the extremum of the entropy, which may be a severe problem. A 
possible explanation to this problem is as follows. The relaxation 
may be incomplete I lLvnden-B elll 1 1 9671) ■ and only the central part 
of system is sufficiently relaxed. Hence the central part is at the 
extremum of the entropy and isothermal, while the outmost part in 
fact can be treated as being isolated from the whole system, and 
thus it is adiabatic. 

Then, if we treat equation ((9) as a differential constraint to the 
above variational calculus, similar to the manipulations of WN87, 
the constrained entropy becomes 



St, r = / F(r,r,,y,y',P,P')dr 



. w dP Gyy 



Anr 2 ' 
dr, 



Gyy' 
2r 



(19) 
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Figure 1. Density profiles of different equations of state. The solid, dashed 
and dotted lines correspond to equations U6t . j2U and 1181 . respectively. 
The long-short dashed line indicates the Einasto profile with a = 0.17. 
Here we take A = —2.27, v = 0.5 and set 4nG = p~2 = r_2 = 1, 
where r_2 is the radius when d In p/d In r = —2. 



where virial theorem has been used because now we have con- 
strained the system to satisfy equation ((9), and 77 is a Lagrangian 
multiplier. Performing a similar variational calculus, we can repro- 
duce the similar result as WN87, 



T) = 



6ivr 2 p 



dr 5dr 7rr 3 



5dr 



(20) 



The problem is that equation 1 ] I St and d20t are both obtained by us- 
ing entropy principle with the constraint of hydrodynamical equi- 
librium, and thus we speculate that the equation of state d20t should 
derive similar density profile as equation (118t . However, KH1 1 pro- 
duces an index-changing polytropic structure, while WN87 exhibits 
a core-halo structure. We provide an explanation to this apparent 
contradiction in Appendix [A] where we also derive an approxi- 
mated solution of equation j20j 



-AP + vP 



A/5 



(21) 



Then equations of state l l 1 !St and l !2 1 1 are found to have similar 
forms, and their solutions can give us similar density profiles which 
are shown in FigureQ] So it is not necessary to introduce the incom- 
plete relaxation here. 

An advantage of equation of state i20\ is that it is at the 
extremum of the entropy, with which we can further study self- 
gravitating systems using statistical-mechanical methods. Then we 
calculate the second order variation of entropy in equation l |19t as 



5 2 St,r 



G 

2r 2 



. G{ ±y ){5 Mf 



WTvr z , .2 67rr 2 p, m 2\. 



(22) 



where the coefficient of (SM) 2 can be positive because A < 0. 
So according to the calculus of variations, our equation (12 1 1 does 
not correspond to the global maximum of the constrained entropy 
of equation H9\ . This result is consistent with one of our ear- 
lier works, in which we find that the equilibrium st ates of self- 
gravi tating systems are the saddle-point entropy states dHe & Kand 
I2OIII) . 
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Figure 2. The surface brightness profile (projected density profile) of 
the elliptical g alaxy. The dashed cur ve comes from the measurement 
of NGC 3379 dCapaccioli et alj Il990l) . and the solid line is calculated 
from equation {21} for a specific choice of the values of A and v. See 
IMo. van den Bosch & W hite (2010) for details about the calculation of p,g. 



2.3 Comparison with observations 

Coincidental ly, equation OH is very similar to the result of a model 
proposed by Ijaffel dl987l) and iHiorth & Madsenl dl993l) . who as- 
sumed that 



N(e) 



for e > 
N for e -> 0" 



(23) 



where N(e) is the differential energy distribution, then the density 
profile is described by an isothermal co re and a poly tropic halo with 
index 5/4. According to the result of Ijaffd dl987l) . this model in 
fact uses an approximation that $ ~ —GM/r at large radii, while 
equation J 1 8fc means that the outmost part is isolated from the sys- 
tem, i.e. $ ~ at large radii, so equation d2H is more accurate tha n 
equation (18). According to the result of lHiorth & Madsenl d 19931) . 
their model can fit the surface photometry along the major axis of 
NGC 3379 very well, and we also present one solution of equa- 
tion J20t in Figure [2] which is compared with the density profile 
of NGC 3379. Besides, as addressed previously, equation J2U can 
be approximated as the trun cated isothermal sphere, which is just 
the King model <Kindll96fj) . then our results are also applicable to 
globular clusters. Hence it is interesting to see that equation J2U 
accounts for observations much better than for the simulation re- 
sults concerning dark matter halos, especially in the center part of 
the system. 



3 FORM OF THE ENTROPY 

In the above section, although we get a density profile with finite 
mass using statistical-mechanical methods, the extremum of en- 
tropy is not a global maximum, which is not consistent with the 
usual principle of maximum entropy. We consider that it may be 
still necessary to obtain equation POt by maximizing some specific 
entropy form, which is constructed from the distribution function /. 

From observa t ions and numeric al simulations (e.g. 

iNavarro et alj 1201 dt iGuimaraes & Sodril 1201 ll) . we know that 
the velocity distributions of stellar systems and dark matter halos 
are anisotropic, while BG entropy maximization with mass and 
energy constraints gives iso t ropic distribution func tion. Some 
studies dBertin & Trenti 20oH iTrenti and Be rtin 2005) considered 



the derivation of anisotropic distribution functions from entropy 
principle. They calculate the extremum of the BG entropy by 
imposing the constraints of mass, energy and a third quantity Q, 
which is related to the single particle's angular momentum and is 
introduced to consider the incomplete violent relaxation. But here 
we do not consider the anisotropy and only impose energy and 
mass constraints on the entropy maximization, so we will derive a 
density profile from an isotropic distribution function. Then notice 
that through the definition of the effective pressure we can always 
get a relation among /?, a 2 and p (see KH1 1), but this relation will 
not be discussed again here. 

It is interesting to notice that the extremum of BG entropy 
leads to an isothermal solution, while the extr emum of Tsall is en- 
tropy gives rise to a polytropic solution (e.g. lHansenll2004T) . and 
equation J2U is just a sum of isothermal and polytropic solutions; 
while BG entropy describes short-range interacting systems and 
Tsallis entropy may describe long-range interacting systems. These 
suggest that the proper entropy for self-gravitating systems, which 
is a functional of /, may be a hybrid of BG entropy and Tsallis 
entropy. Notice that although collisionless gravitating systems are 
described by CBE, we can never neglect collisions at least in the 
central regions (r ~ 0) of the real systems. For these reasons, we 
consider an entropy form as 



S[f]=S[fi,f2]= /(-/iln/i- 



/I ~ h 
(7-1 



)dr, 



(24) 



where /i = ujf, / 2 = (1 — u)f, < ui < 1 and uj is a param- 
eter. Notice that equation MA\ is consistent with equation (2). The 
constraints are: 



/dr = M, 



fedr 



(25) 



where e = v 2 /2 + $. The co nstraint of hydro dynamics equilibrium 
can be realized by / = /(e) dChavani sl2008l) . Then the constrained 
entropy St will be: 



St= / [-/iln/i- 



ft - h 



+a(/i + / 2 )+A(/ 1 +/ 2 )e]dr, (26) 



q-1 

where a is also a Lagrangian multiplier. Evidently, from SS t — 
we get: 

u]n(uf) + {l-u) 9 ^~ U ^" -Ae' = 0, (27) 
q - 1 

where e' = v 2 /2 + $+[(§- u)a - U)]/\(q - 1) = v 2 /2 + 
and we will leave out the prime in the following text. The stabilit y 
of the system requires that /'(e) < dBinnev & Tremaina 2008), 
so A < 0. If w = 1, we can get the result of the BG entropy, so 
p — — \P; if uj approaches 0, from equation ((7} we have 



p = 47r( 



2 (— e) <!-i dv 



(28) 



c„(-$) n = -P~. 



where 



2 n+2 7r(- 



<Z)A, 



7 



3(2 - 



sin f?cos 9- 1 6d6, 



1 - 1 l- 



(29) 



in which 7 is the polytropic index. The upper limit of the integral 
in equation d28t is \] — 2$, because we will ensure that e < for 
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general q. In fact, only we take this upper limit, can we get the last 
step of equation J28b . so we think that this upper limit is responsible 
for finiteness of the mass of the self-gravitating system. Against our 
expectation, we cannot get equation J2U from equation ( 127b when 
q — 7/5 (and hence 7 = 5/4). 

Then we consider another case, that is, the entropy form is the 
same as equation ( 124b . but we assume that fa , fa will be treated 
as independent variables so that 5 St — should hold for fa and 
fa, respectively. With the following two assumptions that Pi = 
P2 = P and p = iop\ + (1 — uj)p2 (P is determined by / through 
equations and ([8}), we derive the result of the extremum of the 
entropy: 

1 X 

pi = —XPi, P2 = -Po , 

K 

so p = -ljXP + (1 - cj)P 1/7 /k = -X'P + u'P lh , which 
must be equivalent to equation ( 121b when q = 7/5. We make these 
assumptions mainly due to the following two reasons: (1) the con- 
dition Pi = P2 = P can ensure E^i = E^2 = Pfc through 
equation d 1 3b . and thus fa, fa and / all satisfy the constraint of the 
energy, i.e. E\ = E2 — E; and (2) p = upi + (1 — uj)p2 indi- 
cates that the spatial density may be a combination of the two terms 
corresponding to fa and fa, respectively. 

For this case the second order variation of the constrained en- 
tropy of equation j26) is 

52 St = \\ {-J^ h) * ~ <ih q - 2 {&f2) 2 )&T, (30) 

which implies that fa increases the BG entropy, fa increases the 
Tsallis entropy, and hence S is globally maximized. 



4 CONCLUSION 

The final statistical equilibrium states of isolated systems can be 
obtained by calculating the extremum of the entropy with the con- 
straints of fixed mass and energy. Based on the dynamical similar- 
ities between self-gravitating systems and fluids, we further con- 
firm the validity of fluid-like entropy for self-gravitating systems. 
Also, we improve the variational calculation of our previous work 
(Kang & He 201 1), and propose an index- varying polytropic struc- 
ture, which has finite mass and can be obtained under the assump- 
tion of incomplete relaxation. If we treat the equation of dynamical 
equilibrium as a dif ferential constraint to the entropy, we derive the 
similar equation as IWhite & Naravanl dl987h . Hence, we conclude 
that with statistical-mechanical methods we can obtain a density 
profile with finite mass without introducing the incomplete relax- 
ation. It is interesting to notice that our results may be more con- 
sistent with observations than numerical simulations. Insisting on 
the validity of the principle of maximum entropy, and based on the 
new equation of state derived in this paper, we suggest a specific 
entropy form, which is just a hybrid of Boltzmann-Gibbs entropy 
and Tsallis entropy. We argue that this hybrid entropy may com- 
pletely describe self-gravitating systems with both short-range and 
long-range interactions. 
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APPENDIX A: 

Here we present our understandings of the solution of equa- 
tion l[20]l, which is different from WN87. WN87 in their figure 2 
showed the numerical solution of equation ( 120) , from which we 
speculate that their solution has no essential difference from the 
cored isothermal sphere (p = — AP), which is an incomplete solu- 
tion and has infinite mass. But we consider that isothermal sphere is 
not the only solution of equation ( 120b because of two mathematical 
reasons: 

(i) We set m(r) = 3r~ 3 r 2 p/ Pdr = n(r)/2nr 3 . If 
m(r) = const, we must have p/P = const. Then from equa- 
tion 120b we know p — — \P and m = — A. However, there is no 
physical requirement that m(r) = const, i.e. m(r) is commonly 
not a constant, so more general solutions must exist. 

(ii) We can rewrite equation ( |20b as one differential equation by 
eliminating n(r) 

5^P'(r)r + 15^+3A = 12£ J (Al) 

where P'(r) is the derivative of P with respect to r. Because 
P = P(p, r), we can treat r and P'(r) as functions of p and P. 
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Figure Al. The relative error between m and p/P as a function of r. 

Then equation JAU is a second-order differential equation of p with 
respect to P, so there must be two integral constants in the solution 
p — p(P). The only physically reasonable boundary condition that 
we used in this paper is P = when p = 0, which was already 
used by WN87, and thus there is only one integral constant left. 
However, WN87 also set P| r= o = p\r=o = 1, which eliminates 
the only remaining integral constant and is not proved to be reason- 
able. This is the reason why WN87 "comes to an important and sur- 
prising point" and "contains no model with a concentration larger 
than 3.35". WN87 considered that another shape parameter should 
be the minimal radius they introduced, but this quantity may never 
be treated as an integral constant. Hence the solution p = p(P) 
must have two shape parameters, including a Lagrangian multiplier 
A and an unknown integral constant. 

Because m(r) is an arbitrary function of r, the general exact so- 
lutions are very difficult to obtain, and thus we consider to use 
the iterative method to look for the solution. We first assume that 
m(r) ~ p(r)/P(r), since (1) when r — > and r — > oo, then 
m(r) — > p/P according to the L'Hopital's rule; and (2) m(r) — 
p(r)/P(r) can exactly set up only when m = —A = p/P; so 
when m(r) is not a constant, p(r)/P(r) can be a more accurate 
replacement of m(r). Then we can approximately rewrite the equa- 
tion i20\ as 

2- -£<*-*£>■ 

whose solution is 

p = -\P + vP 4/5 , (A3) 

where v just is an integral constant. For the general value of v, if 
one wants to get a more exact solution, he can first denote dA3t as 

p {0) =p (o)(A,^p),then 

(i) substitute p' ' into m(r), which can be rewritten as m(r) = 
m(p, P) = m(A, v, P), then from equation HQ) we can get a new 
density p (1) = p' 1 ' (A, v, P); 

(ii) repeat the step of (i), until the maximum of |p'" +1 '(P) — 
p^ n '{P) \ becomes smaller than e, which is the resolution we set. 

Notice that every time we iterate the function p = p(A, v, P), not 
the value of p at a point. Then p' n ' will be the solution of equa- 
tion J20t with the resolution e. 

However, equation dA3t has the main properties of equa- 
tion ( |20| (. In fact, m ~ p/P is really a good approximation except 
for large value of v, and hence the typical error of m is small for 
small value of v. We define the relative error between m and p/P 



as R = \(rn — p/P)/m\ = [1 — p/Prn\, which is a function of 
r. We show R vs. r in Figure IaTI in which all the data are taken 
from the results corresponding to equation dA3l , i.e. equation J2U . 
in Figure[T] From Figure[TJwe know that when r/r_2 < 0.2, equa- 
tion ( IA3b can be very well fitted by p/P = —A, which almost 
ensures R = 0, so we only focus on the external part. Notice that 
when r > r_2, both p and P are very small, so even a small devia- 
tion from p or P may cause a high value of R. Figure lATl indicates 
that with a chosen value of v, the approximation we adopted is 
very good. If we make more accurate approximation of m(r), we 
can obtain more accurate solution of equation i20\ , but the shape 
of density profile will not be changed essentially. 



